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Abstract 



^T) . We present a number of analytical results which should guide the interpretation 

^^ i of lattice data in theories with an infra-red fixed point (IRFP) deformed by a mass 

^O [ term 6C = —mqq. From renormalization group (RG) arguments we obtain the 

C^ I leading scaling exponent, F ~ m''^, for all decay constants of the lowest lying 

states other than the ones affected by the chiral anomaly and the tensor ones. These 
scaling relations provide a clear cut way to distinguish a theory with an IRFP from 
^ ■ a confining theory with heavy fermions. Moreover, we present a derivation relating 

the scaling of (qq) ~ m''''' to the scaling of the density of eigenvalues of the massless 
Dirac operator p(A) ~ A''"''. RG arguments yield r]gg = (3— 7*)/(H-7*) as a function 
of the mass anomalous dimension 7* at the IRFP. The arguments can be generalized 
to other condensates such as (G^) ~ m^'^^'^'^*' . We describe a heuristic derivation 
of the result on the condensates, which provides interesting connections between 
different approaches. Our results are compared with existing data from numerical 
studies of SU(2) with two adjoint Dirac fermions. 
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1 Introduction 

There are numerous examples of two-dimensional field theories that are invariant under 
the full conformal group. In four dimensions, the beta function of A/" = 4 super Yang-Mills 
is known to vanish to all orders in perturbation theory, for any value of the coupling, so 
that the theory is scale invariant. Other theories have isolated zeroes of the beta function 
that correspond to fixed points of the renormalization group (RG) flow. For instance, 
the gauge coupling g in QCD flows to zero as the energy scale is increased, leading to 
the well-known phenomenon of asymptotic freedom; in this case (7 = is commonly 
called an UV fixed point. On the other hand, if a theory has an IR fixed point (IRFP), 
the couplings will flow to such a fixed point at large distances, and the theory becomes 
scale invariant in the large-distance regime. Theories with an IRFP do not break chiral 
symmetry spontaneously, and are said to lie in the conformal window. 

Supersymmetric examples of theories within the conformal window have been stud- 
ied in detail - see e.g. Ref. [1] for a review. Recently there has been a lot of interest 
in identifying non-supersymmetric gauge theories with an IRFP. The main motivation, 
besides intrinsic interest, comes from the fact that theories near the conformal window 
correspond to the class of theories underlying walking technicolor [21 El IH El [6] , which is 
the phenomenologically most viable offspring of technicolor theories O [H |9l [10]. A chi- 
rally broken theory near the edge of the conformal window is supposedly identified by an 
enhancement of the ratio {qq)/f^ with respect to a QCD-like theory [10]. Unfortunately 
this quantity does not display a simple known parametric behaviour. Another strategy, 
which is adopted in this paper, is to first identify theories within the conformal window, 
and then approach the boundary of the window using the available information on the 
color- flavor phase diagram [21 [3l [6l [11] . 

The identification of conformal theories using numerical simulations is a difficult task, 
since the only observable quantities would be the power-law scaling of correlators at large 
distances. However actual lattice simulations are performed in a finite volume, and with 
a non- vanishing fermion mass; both the mass and the finite size of the system are relevant 
operators at large distances and drive the theory away from conformal behaviour. Turning 
a technical limitation into a tool, it has become a standard strategy to consider conformal 
gauge theories (CGT) candidates deformed by a mass term, and to identify them from 
the study of their hadronic observables. Thus, if there exists an IRFP, the lattice results 
should be described by a mass-deformed conformal gauge theory (mCGT), obtained by 
adding a bare mass to the original lagrangian 

6C = —mqq. (1) 

As a consequence of the deformation, these theories are expected to develop a mass gap 
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and a fermion condensate and thus give rise to asymptotic states and related observables, 
which scale to zero as the massless limit is approached. For any observable O the leading 
exponent rja of the mass deformation is defined from its scaling as rra — t- 0: 

O ~ m'^° + higher order in m + terms analytic in m . (2) 

These critical exponents can be measured on the lattice and it is the aim of this work to 
provide predictions for them that can be tested numerically. 

The paper is organized as follows. In section 11.11 we set the framework by discussing 
some characteristics of theories inside the conformal window. In section [2] we discuss 
general aspects of IRFPs, and introduce the standard tools for analyzing the behaviour 
of field correlators near a fixed point of the RG flow. Thereby we obtain the hyperscaling 
relations that are usually derived in the context of critical phenomena [12], and we study 
the information that they yield in the framework of mCGT. 

Section |3] is devoted to the study of the chiral condensate in mCGTs. First we review 
the relation between the scaling of the chiral condensate with the fermion mass, and the 
density of eigenvalues of the massless Dirac operator in the infinite- volume limit. As 
stated above, the chiral condensate must vanish as the fermion mass is taken to zero at 
a rate that is dictated by a critical exponent r]qq. The non-analytic dependence of the 
fermion condensate on the fermion mass is directly related to the scaling exponent for 
the eigenvalue density of the massless Dirac operator. As pointed out in Ref. [13], the 
exponents turn out to be the same: 

{qq) ~ m"?' => p{X) ~ A^'« . (3) 

The scaling exponent rjqq is determined as a function of 7^,, the anomalous dimension of 
the mass at the IRFP. The RG analysis, which applies to all condensates, yields 

(3 - 7.) .4. 

^'^ (1 + 7*) ^' 

We then present the determination of this coefficient from a heuristic calculation, which 
provides some physical insight in the dynamics of mCGT. The limitations of such a 
heuristic approach are highlighted, and the interpretation of IR and UV cutoffs is clarified. 
We conclude this section by analyzing current lattice data for the eigenvalue distribution 
in an SU(2) gauge theory with two flavours in the adjoint representation. 

In section H] we explore the consequences of hyperscaling for the decay constants of 
the hadronic states. Our results, summarized in Tab. [H can schematically written as, 

G^m^, {0\O{0)\H{p)) = G , (5) 



for operators with scaling dimension A. Further informations are obtained by combining 
these results with the chiral Ward identities in section I4.2( these scaling predictions for 
the decay constants are then compared with recent results from numerical simulations of 
potential mCGT on the lattice. Finally we discuss the implication of the scaling of the 
decay constants for the width of the hadronic states, and compare the scaling of the decay 
constants in a mCGT to the one of heavy quarkonia states in a chirally broken theory 
like QCD. 

1.1 Conformal window - discussion and results 

It is well known that SU(A^) gauge theories with Uf fermions are asymptotically free as 
long as Uf does not exceed an upper hmit that depends on the number of colours Nc 
and the fermion representation R. At small distances the gauge coupling decreases log- 
arithmically, and the dynamics is successfully described by perturbation theory. In the 
SU(3) gauge theory minimally coupled to Uf = 2 light flavors in the fundamental repre- 
sentation, the coupling increases at large distances, and the theory undergoes confinement 
and spontaneous chiral symmetry breaking, exhibiting a spectrum of bound states. In 
the massless limit, the spectrum includes three massless Goldstone bosons, known as 
(vr", TT"*", 7r~), reflecting the spontaneous breaking of chiral symmetry. As a consequence, 
there is a gap in the spectrum between the pions and the rest of the states whose masses 
are parametrically of the order of some hadronic scale A ~ Aqcd? and remain flnite in the 
chiral limit. At low energies compared to A the dynamics are successfully described by 
an effective theory of self-interacting pions, known as chiral perturbation theory. A small 
non- vanishing mass can easily be incorporated as a perturbation of the massless theory. 
As the number of light fermions is increased, before asymptotic freedom is lost, the 
theory may develop an infrared flxed point (IRFP) due to the effect of the fermions 
on the running of the coupling. We shall denote by nj^c the number of fermions above 
which the theory exhibits an IRFP. In this case the theory becomes scale-invariant at 
large distances, while the short-distance behaviour is still the one dictated by asymptotic 
freedom. As a consequence of the scale invariance at large distances, the theory cannot be 
in a conflning phase and chiral symmetry remains unbroken. The long-distance dynamics 
is governed by the critical exponents of the IRFP, which determine the scaling laws in the 
vicinity of the flxed point. The Banks-Zaks theories [H], where Nc and n/ are arranged 
such that the critical coupling g* '^ 1, provide one working example of a theory within 
the conformal window. Early studies of near conformal and IRFP theories were based 
on approximate solutions of the Schwinger-Dyson equations [T^l [ISj ; these analyses were 
extended to higher representations in Ref . [11] . Unfortunately it is very difficult to control 



the systematic errors due to the truncation of the IPI vertices appearing in the Schwinger- 
Dyson equations. Moreover Schwinger-Dyson equations predict the anomalous dimension 
of the mass to be around one, whereas unitarity constraints on the conformal group [17], 
in principle, allow for 7* < 2. 

Recent results have appeared recently, that address this problem either from an RG 
point of view [HI UHl |20], or from a gauge/string duality perspective [211 1221 [23]. We 
defer the investigation of the connections between our results and these other approaches 
for further studies. 

Recent numerical simulations of gauge theories on the lattice have triggered a renewed 
interest in those theories and in turn in technicolor models. Algorithmic progresses have 
made lattice simulations with light dynamical fermions accessible on current hardware [2H 
[25l [26] . This opens the possibility to obtain first principles results for technicolor, and 
several preliminary investigations have appeared [271 [2H1 [29l [301 [311 [321 [331 [311 [33 [33 
[371[3Hl[3Sl[ini[I31[lIl[121[131[Sl[lS]. It is important to bare in mind that recent lattice 
results for theories that may lie inside the conformal window are plagued by systematic 
errors, and their interpretation still needs to be clarified. A recent discussion of the lattice 
artefacts in simulations of theories with a potential IRFP can be found in Refs. [lH [46] . 
For these theories, unlike in QCD, there are no experimental results to guide the lattice 
simulations. 

Therefore it is crucial to develop analytical results in order to guide the lattice studies, 
and help in analyzing their outcome. A wider range of analytical predictions, together 
with more extensive simulations, will help in finding robust evidence for the existence of 
IRFPs. 

2 Infrared fixed points 

Let us henceforth consider theories inside the conformal window, i.e. gauge theories 
minimally coupled to a number Uf of Dirac fermions, with the number of flavors and their 
representation adjusted so that the theories are scale-invariant at large distances when 
the fermions are massless. In general, fixed points of RG flows are identified by the zeroes 
of the P functions that describe the evolution of dimensionless couplings. The typical 
evolution of a running coupling is sketched in Fig. [H The running coupling flows to a 
constant value at small energies, which corresponds to a zero of the beta function. The 
value g* of the coupling at the fixed point, and the precise shape of the nonperturbative 
function g{fi) are scheme- dependent. However the existence of the fixed point and the 
critical exponents are universal. 

The fermion mass is a relevant coupling at the IRFP, and drives the theory away 
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Figure 1 : Running of the coupling as a function of the energy scale for a theory with an IRFP. 
At low energies the coupling flows to a fixed-point value g* , while the high energy behaviour 
is the usual one expected for asymptotically free theories. The scale Ku corresponds to the 
energy where the running starts to be dictated by asymptotic freedom. The dashed curve at 
low energies shows the running of the coupling when a fermionic mass term is switched on. 



from it. In a theory with a non- vanishing fermion mass, the fermionic degrees of freedom 
decouple at low energies, and the theory behaves like a pure Yang-Mills theory. The 
running of the gauge coupling for the massive theory is given by the dashed curve at 
small /i in Fig. [1], where the running of the coupling below some scale Air is explicitly 
drawn. Note that in the presence of an IRFP Air goes to zero as the fermion mass 
vanishes. 

The running of the mass is described by its anomalous dimension, which has the 
opposite sign of the anomalous dimension of the renormalized composite operator gg. 



d _ , 
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We have explicitly indicated the scale dependence of the various quantities. In this paper 
we will use the symbol 7 to denote the anomalous dimension of the mass and quark 
condensate: 7 = 7^ = —'Jqq. 

Note that the anomalous dimension away from the fixed point depends on the renor- 
malization scheme. However its value 7^. at the IRFP is a scheme- independent quantity. 
A concise discussion of the scheme- dependent features of IRFPs can be found in Ref. |16]. 

Throughout this paper we will often refer to scaling dimensions of operators, denoted 
by A; thery are obtained as the sum of the naive mass dimension of the operator and the 



anomalous dimension. For example for the operator qq we write: 

^qq = dgg + -fgg = 3 - 'y^ , 7/^ = 1 + 7=^ , (7) 

where we have also introduced the scaling exponent ym, which often appears in what 
follows and is widely used in the RG-literature [12]. Throughout this paper we will use 
these notations interchangeably. 

Scaling laws are derived by assuming that the fermion mass is the only relevant op- 
erator at the IRFP. RG equations will be used below in order to derive the scaling of 
the chiral condensate as a function of the fermion mass. It is therefore worthwhile to 
briefly recall how the scaling relation for the masses in the spectrum is obtained. A recent 
discussion of RG flows in the vicinity of an IRFP can be found in Refs. [101 SZl SB] • 

Let us consider the zero-momentum vacuum correlator of an interpolating field H{x) 
with the quantum numbers of a given state in the spectrum: 



CH{t;g,rh,fi) = j d\ {H{t,x)H{Q)\.^^^ 



where we have indicated explicitly the dependence on the couplings and the scale ii. It 
is useful in this context to introduce a rescaled mass m(/i) = m(/x)//i. For the specific 
case of lattice simulations, the scale is set by the inverse lattice spacing /^ = a~^ . The 
masses of the physical stable states are obtained from the Euclidean time dependence of 
two-point functions. At large Euclidean time t: 

CH{t-g,m,fi)^e-^'"\ (9) 

where Mh is the mass of the lightest state in the channel under examination. We examine 
the consequences of the RG equation for the two-point function. 

In the vicinity of the fixed point, a RG transformation acts on the correlator according 
to: 

/i = V; C^(t;(?,m,/i) = r2^«C^(t;(7',m',/i'), (10) 

where jh is the anomalous dimension of the field H. The flow of the couplings near the 
RG fixed point is power-like: 

g' = W^g , m = W'^rh . (11) 

We shall neglect henceforth the irrelevant coupling g {yg < 0). Multiplying all mass units 
by the factor h we obtain: 

Cuit] m', fi') = h'^'^^Cnith-^] m', fi) , (12) 



where dn is the naive mass dimension of the operator H . Choosing h such that m' = 1, 
the equations above yield: 

CH(t;m,/i) = C^F(tm^/(i+>),/i) , (13) 

where F is some function that, for fixed //, depends on the rescaled variable x = tm^^^^'^'^*'' 
only. The detailed dependence of the prefactor Ch on the parameters of the theory is 
postponed to the next section, where it will play a prominent role. Comparing Eq. (1131) 
with the expected behaviour Eq. (jH]) yields: 

Mh ~ CHjJi.rh^+-i* as m — )■ . (14) 

Note that the scaling of the mass Mh is entirely determined by the anomalous dimension 
7* and does not depend on the specific choice of the interpolating operator H. Eq. (IT^ 
shows that all lowest state masses scale with with same exponent 1/(1 + 7*), while the 
proportionality constant ch depends on the chosen channel. While each individual mass 
in the spectrum vanishes, ratios of masses should remain constant as the chiral limit is 
approached. This scaling is consistent e.g. with the scenarios proposed in Ref. liniEO]- 

In the derivation above we have not considered the effects of a finite decay width. At 
least one channel ought to be stable and therefore not affected by the width. According to 
an inequality by Weingarten [51], valid for np > 2, this should be the mass of the lowest 
pseudoscalar flavour- nonsinglet, which we shall later on denote by Mpa. For all other 
states one might wonder how the width interferes with the derivation above. Could the 
width and the mass conspire to cancel their leading mass scaling behaviour in such a way 
as to invalidate Eq. flT^ ? We would like to bring forward two reasons why this should 
not be the case. First the difference in the large A^c-scaling of mass and width [Th/Mh ~ 
0{1/Nc)) from QCD should hold in mCGT too and serve as a parametric argument against 
such a cancellation. Second we show in appendix O that in the approximation where the 
self-energy is treated as being constant such a cancellation can be excluded. This seems 
intuitively plausible since in Euclidian time the mass and decay width behaviour are 
associated with exponential and oscillatory behaviour respectively. 

On the contrary since mass and width do not seem to interfere in the leading large 
i-behaviour Eq. (1131) suggests that both the mass and the width of the resonance scale 
according to 

M,r~mi/(^+^*), (15) 

We shall revisit the scaling of the width in Sect. HJ after discussing the scaling of the decay 
constants and derive T{A ^ B + C) ~ rn^/i'^+'y*) for a specific decay A ^ B + C. 

The behaviour ( !T^ is markedly different from what is observed in the spectrum of 
theories where chiral symmetry is spontaneously broken, like e.g. in QCD. In the latter 



theories, the Goldstone bosons become massless in the chiral hmit, while the other states 
remain massive, with their masses being of the order of some typical hadronic scale A. For 
theories with an IRFP, all states become massless, presumably at the same rate, which 
prevents a simple description of the nonperturbative low energy hadronic dynamics in 
terms of an effective theory like chiral perturbation theory. 

Let us conclude this section by recalling how the finite-size effects can be analyzed 
using RG equations. We shall discuss explicitly the case of the correlator Ch, including 
the dependence on the size of the system L. We remind the reader that by studying finite 
volume effects, it is implied that the box is larger than the typical scale, L ^ /i~^, and 
therefore does not interfere with characteristic short distance dynamics. The solution of 
the RG equation, including the L-dependence, scales as, 

Cnit- m, L, /i) = h-^^"CH{t- m', L, //') , (16) 

according to a modified version of Eq. (ITOj) . Rescaling the energies by the factor 6, and 
using the power-law scaling of the couplings near the IRFP yields: 

CH{t] m, L, /i) = b-^^'^"+^»^CH{b-H; W^m, b'^L, //) . (17) 

Choosing b such that b^^L = Lq, where Lq is a reference length, yields: 



where we have introduced the scaling variable x = Ly^^m. 

Comparing Eq. (fT8|) with the expected asymptotic behaviour in Eq. IQ we obtain: 

Mh = L-'f{x) , (19) 

where f{x) is some function of the scaling variable x, expected to vanish when x goes to 
zero. In order to recover the correct scaling with m in the thermodynamic limit 

/(x) ~ x^/^'" , as X -> oo . (20) 

As one can see from Eq. ( TT9|) . if the fermion mass is decreased at fixed /x and L, then the 
mass of the states in the spectrum will initially decrease until the Compton wavelength 
of the states is of the order of the linear size of the system. When this happens, the mass 
of the states saturates and scales with the inverse size L~^. Results for MhL computed 
on different volumes should follow a universal curve when studied as a function of the 
scaling variable x. 
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3 Modified Banks-Casher relation 

In this section we relate the scahng exponent of the chiral condensate r]qq to the scahng 
of the eigenvalue density of the massless Dirac operator. We then illustrate how the 
RG equations yield a prediction for the exponent in terms of the anomalous dimension 
7* introduced in Eq. ([6]). These results follow readily from the RG scaling of the free 
energy and the field correlators in the vicinity of fixed point, and were already presented 
in Ref. [I3]. Here we discuss in detail the derivation of these results in the context of 
a mCGT, generalizing to other condensates such as the gluon condensate, and compare 
them to a more heuristic derivation. 

3.1 Eigenvalue density p{X) and the scaling exponent rjqq 

It is useful to recall the basic steps in the derivation of the Banks-Casher formula, in order 
to highlight the order in which the limits are taken, the divergences that may appear, and 
to identify the differences from the case of a conformal theory. 

We closely follow the discussion in Ref. [52] and extend it at appropriate places to 
mCGT. The fermion propagator can be written as: 

(,W5(,)) = 5;^^^Wii!M. (21) 

^^ m — iXn 

n 

where the eigenmodes of the massless Euclidean operator D = "y^D^ have been introduced: 

Dun{x) = A„M„(x) . (22) 

Since the eigenfunctions occur in pairs with opposite eigenvalues, the chiral condensate 
in a finite volume V is given by: 

{Q(1)v = — dx (g(x)g(x)) = — — V . (23) 

VJ V ^^m^ + Xl 

Taking the infinite volume limit at fixed mass, the sum over positive eigenvalues can be 
replaced by: 

P(A) 



oo 



{qq) = lim {qq)v = -2m / dX ;y , (24) 

where p(A) denotes the number density of eigenvalues per unit volume. Eq. (1241) is purely 
formal at this stage in the sense that a UV-regularization is needed on both sides. In four 
dimension the divergences are logarithmic arrd quadratic respecfvelJH The divergences 



^Note that if the regulated theory breaks chiral symmetry explicitly, as is the case with lattice Wilson 
fcrmions, then a cubic divergence appears that survives in the chiral limit [53j . 
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can be isolated via a twice-subtracted spectral representation: 



" PW o,„5 r ''>• PW . ... ^3 



m = -2m j^ dX-,^^ - 2m' j^ -^^j^ + T.m + -,,m' . (25) 

The subtraction constants 71 and 72 contain the UV- divergences. Their respective be- 
haviours are 71 ~ A^y, and 72 ~ log [A^v], and their actual values depend on two physical 
renormalization conditions used to define the finite condensate on the LHS of Eq. (124|) . 
We shall investigate the limiting behaviour when m — ?> 0. The second integral and the 
subtraction terms in Eq. f l25p vanish in the chiral limit {m — > 0). Therefore only the 
first integral, sensitive to the IR region, can result in a non-analytic term and has to be 
investigated further. A simple change of variable yields: 

{qq) = -2 / dx ^^ + A{m) , (26) 

Jo 1 + x^ 

where ^(rra) stands for an analytic function of m. From Eq. (l26|) . following the same 
arguments used in QCD, one can readily obtain: 

{qq) '"~° m''^-' ^ p(A) ^~° A^^^' . (27) 

This in turn implies: 

^7^^00-11^ = , '7wlmCGT > , (28) 

since in QCD the condensate remains finite in the chiral limit, while it vanishes in mCGT. 
Let us derive the same scaling coefficient 77^5 (jl]) from a RG analysis. The starting point 
is the two-point function Cgg{t] m, /x), as in Eq. ([8]), where the hadronic field H = qq, and 
the explicit dependence on the coupling g is suppressed. The solution of the RG equations 
for this specific case is: 

Cgg{t; m, /i) = r2^-C,-,(tri; 6^™m, fi) . (29) 

Imposing again i^^rh = 1, finally leads to: 

Cgg{t; m, /i) = m^^Cgg{tm^/y-; 1, /i) . (30) 

Inserting a complete set of states the exponential decrease of any state other than the 
vacuum for large t results in: 

C,g{t;m,fi)'^°°m'^^^ , (31) 

whence the scaling exponent (!27|) follows: 

--^-^- '-) 
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The eigenvalue density then scales as: 

p(A) ~ Xi^-^*yi^+^*) ; (33) 

this result generalizes the Banks-Casher relation for QCD ^^■. 

m L=o ^ => p(0) = -TT (qq) L^o (34) 



to mCGT. It is interesting to remark that Refs. [13] and [551 HS] state different predictions 
for the scaling exponent. Our determination of this critical exponent agrees with Ref. [T3] . 
Surely this derivation generalizes to any other operator, for example the gluon con- 
densate for which one gets: 

Vg^ = = ^ • (35) 

Um 1+7* 

The scaling dimension of the gluon condensate is four since it appears in the Lagrangian 
density of a four dimensional scale invariant theory. 

3.2 Alternative and heuristic derivation of rjqq 

Let us now present an alternative derivation of the scaling exponents rjqq and 77^2 in 
Eqs. (132 p . (I35p . which is of a heuristic nature but might provide some physical insight. The 
discussion for (qq), which we shall adopt here before generalizing it to (G^) closely follows 
Ref. |56)i rl. In this work we refine the discussion and interpretation of IR and UV-terms 
by making use of the scaling of the hadronic masses in Eq. ( fT4l) and the interpretation of 
subtraction terms in Eq. ( 12^ . 

In a low energy effective theory describing the dynamics of the operator qq □, the mass 
deformation in Eq. ([1]) corresponds to a tadpole term and demands a reminimization of 
the potential to find the stable vacuum. The potential for qq is not known but the scaling 
of the two-point function is governed by the anomalous dimension. It has been proposed 
in Ref. [SS] to mimic the continuous spectrum of such an operator by introducing a tower 
of scalar fields with suitably adjusted masses and couplings: 






qq{x) r^ 2^ fnipnix) ; (v'nIgglO) ~ /„, , \ ni2 _ ,2 ' f^^^) 




S' {Ml 



^aa 



^Tlie computation in Ref. [56] differs by in an additional term 5C ^ {qq)^ which is not relevant here. 

■^The calculation is similar to an analysis of a scale invariant theory with a scalar operator and tadpole 
term Ref. |57j in the context of the unparticle scenario, where 2 < 7* < 1 (A;^ = 8 — 7*) was assumed 
and made it necessary to introduce (various) IR regularizations. 

^Wc refrain to change to a notation qq — > Ojj since we are not interested in parametrizing an effective 
theory for Ojj as in Ref. |56j . 
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where the quantity S describes the mass spacing between the v9„-modes. The decom- 
position (136!) reproduces the two-point function of a conformal theory in Minkowski 
space in the hmit 5 — )■ [58], up to potential subtraction ambiguities. Note that in 
Eq. (l36l) we have not tried to keep track of the overall mass dimension and normaliza- 
tion since they are irrelevant for scaling properties. The potential part of the Lagrangian 
C = —m'^^fnipn — '^/'^Ylin^n'^'n then leads to the equation of motion for </?„ of the 
form: 

mfn + Ml^n = ^ (v9„) = -mfn/Ml , (37) 

with solution as indicated on the right. Thus leading to a VEV, 



n " 

where IR- and UV-cutoffs, to be discussed below, were introduced. We have taken into 
account the running of the gauge coupling as indicated in Fig. [H though showing a 
somewhat cavalier attitude towards the treatment of the transition region to be justified 
later on. The non- analytic part in m, if present, is hidden in the IR-cutoff. It seems 
natural that the latter is governed by the typical hadronic mass scale, i.e. Air ^ cMh, 
where c is a constant irrelevant to our investigations. The integral can be computed and 
yields: 

{qq) ~ -m (M^)^'-'"' + miA'X"^-' A^^ . (39) 

Thus using the scaling of the hadronic masses ( !T4|) and ([7]), Eq. ( 139|) becomes 

{qq) ~ m^ + A{m) => rjqg = -—^ , (40) 

1 + 7* 

where as previously A{m) ~ 0{m) denotes an analytic function in m. We have therefore 
derived the exponent riqq in Eq. ([3]). The UV-divergent term in Eq. (!39|) corresponds 
to the quadratic divergence discussed in the previous section and is irrelevant for the 
non- analytic part and the scaling exponent rjgg. The non-appearance of the logarithmic 
divergence might be related to the fact that we do not consider the back reaction of the 
mass perturbation on the spectrum such as taking into consideration power correction in 
m in the couplings /„. 

Surely this procedure generalizes to any gauge invariant term in the Largangian 6C = 
^i-Ao)/ymQ jj^ which case the condensate (HOl) assumes the form: 

(C) liR ~ m^ (M^)^o~2 _ ^1^ ^ (^4;^) 
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in accordance with Eq. (!32l) which is absolutely general. The subscript IR indicates that 
UV-terms have been omitted. 

We consider it worthwhile to discuss the term 5C, ~ G^, resulting in the gluon con- 
densate. From Eq. (HTl) . paying attention to the UV-terms in addition one gets: 

and with (IT^ and Aq2 = 4 one gets : 

(G2)~m^^^SF^+^(mO) => r^G^ = ^^^^' ~ ^^ = -^- , (43) 

1/m 1+7* 

the same scaling as in Eq. (j35])- The A{'wP) refers to the term proportional to A^y. It 
originates from the region of asymptotic freedom and is interpreted as a mixing with the 
identity. Such contributions, sometimes called renormalons, have hitherto prevented an 
extraction or a proper definition of a gluon condensate in lattice QCD. 

3.3 Critical discussion 

Although the derivation above reproduces the correct result there remain some points that 
deserve further clarification. For this discussion we shall first think of the theory when 
m = 0. Assuming that this effective field theory approach, the deconstruction (15^ . can 
be extended to higher orders one would need to introduce higher order terms, starting 
from cubic ones, into the effective Lagrangian in order to reproduce higher correlation 
functions. There are two important effects, due to higher order terms. 

First, are they going to modify the leading order extraction of the {qq) directly? 
The answer appears to be no, since the initial potential does not know about the small 
parameter m and perturbing the system by the term in Eq. ([1]) does not lead to major 
correction to any order but in the linear one. That this is self-consistent can be seen 
more explicitly by plugging in the VEV (</?„) ~ ?7i as in Eq. f l37p into a fictitious higher 
order term. Second, does the quadratic order need to be modified? Since the higher order 
terms are going to modify the two-point function, the answer appears to be yes. As long 
as these modifications are of the form M^| higher order ~ n'^ with a > 1, they will give rise 
to subleading effects in m. Although this appears likely we have not tried to resolve this 
issue in this paper but obviously this question deserves further study |j. 



■'One might also be eoncerned that the mixing of different (^„ modes might reshuffle hierarchies. That 
this is not the case for the analytic part follows from the fact that the modes with low n are responsible 
for the non-analytic part and that the higher modes are numerically supressed w.r.t to the lower modes 
by (^„) = (^i)/ni+^- m- 
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3.4 Lattice data 

Note that the anomalous dimension 7^, is related to the running of the fermion mass, and 
has been computed by Schrodinger Functional methods in Ref. [16]. On the other hand, 
as discussed above, the exponent riqq characterizes the behaviour of the eigenvalue density 
around zero, so that in principle it can be extracted from the eigenvalue density. 

First results for the eigenvalue spectrum of the Dirac operator have been presented 
in Refs. [32l 147] for an SU(2) gauge theory with two Dirac fermions in the adjoint repre- 
sentation. An extensive study of the 200 lowest eigenvalues is available only for the 16^ 
lattice studied in the references above. As argued in Ref. [59], the mode number of the 
massive Dirac operator: 



u{M,m) 



+A 



dXp{X) 



(44) 



where A = a/M^" 



m^ 



carries the same information as the density itself. The mode 



number can be renormalized and yields a RG invariant, universal quantity that describes 
the physical properties of the Dirac spectrum independent of the regularization used. 

The mode number can easily be computed from the available eigenvalue distributions. 
Results obtained at /3 = 2.25 on a 32 x 16^ lattice are reported in Fig. |2] for different 
values of the fermion ation of the mode 
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Figure 2: The mode number ^{M, m) for the SU(2) gauge theory with two fermions in the 
adjoint representation. The data show the dependence of the mode number on the scale A for 
several values of the quark mass. 

number to the thermodynamic and massless limit: 



lim lim u{M,m) 

m— >0 y— >oo 



(45) 
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Using the current data at a single value of the volume, where a reasonable amount of data 
is available Ref. [17], the first extrapolation cannot be performed. We defer this analysis 
for further studies as larger lattices become available, and concentrate instead on the 
extrapolation to the chiral limit. The data in Fig. [2] show that there is a dependence of 
the mode number on the PCAC mass. The data at the two lightest masses are compatible 
within the statistical errors, and we shall take the data at the lightest mass for our analysis. 
In the chiral limit, the mode number is expected to scale as: 



u{K,^)=C{K-g) 



'799+1 



(46) 



where the possibility of a non- vanishing spectral gap g, as suggested for instance in 
Ref. [SD] , is taken into account in the functional form of the fitting function, we therefore 
end up with fitting the data to three parameters, namely C, (?, and the exponent. The 
data are consistent with such a power-law behaviour, however the value of the exponent 
depends critically on the range used for the fit. The plot in the left panel of Fig. |3] shows 
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Figure 3: Fit of the mode number v{M^m) to a power law behaviour. The two plots represent 
two fits to the same data over two different ranges. 

the result of a fit that has a fit range A G [0.002, 0.014], in lattice units. The result for the 
critical exponent is 1 + ri^q = 1.65(4), which can be translated into a fitted value for the 
mass anomalous dimension, yielding 7^, = 1.42(5). Note that this value is quite different 
from the one found in previous studies [171 H5] . However, if the fit range is reduced to 
A G [0.003,0.008], then the fitted exponent is 7* = 1.1(3). The latter value is still larger 
than the one obtained form the Schrodinger functional studies, it clearly shows that the 
fitted exponent depends critically on the fit range. More extensive data on the eigenvalue 
distributions are needed in order to be able to extract the critical exponent in a reliable 
manner. 
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Note that the error induced by varying the fitting range turns out to be larger than 
the statistical error, and that this result is obtained at one value of the lattice volume, 
and could be affected by finite size effects. A more comprehensive analysis of the volume 
dependence of the eigenvalue distribution is needed. More extensive lattice data should 
become available in the near future. 

4 Decay constants 

This section explores the information that can be gathered from the RG scaling of matrix 
elements of given operators. In subsection 14.11 we show that the scaling of decay con- 
stants of the lowest lying states directly follow from its anomalous dimension through the 
Callan-Symanzik equations. In subsection 14.21 we deduce consequences from spectral rep- 
resentations of the Ward identities and low energy theorems for pseudoscalar and scalar 
states evaluated at zero momentum. In subsection 14.31 we compare our theoretical predic- 
tions with recent lattice data. Finally miscellaneous matters of interest are presented in 
section 14.41 In appendix |X] the pseudoscalar WI at large momentum transfer is used to 
deduce further information. 

4.1 Hyperscaling and decay constants 

Let us consider an operator O with scaling dimension Aq and quantum numbers such 
that it couples to a state \H{p)) with strength Gh for scalar operators and Fh for vector 
operators. We shall choose Gh to exemplify the equations below: 

{0\OmH{p)) = GH, Ao = do + io. (47) 

Information on the lowest lying state can be gained from the large time behaviour of the 
Euclidian two-point function Co(t; g,m, fi) defined in Eq. ([S]): 



i— ^oo 



_,,^,{0\OmH{p)){H{p)\O{0)\0) _ ,,^, \Gh 



2 



C„(, ,, ™, ,) -^ e-".' -— ^ e-"».^^ (48) 

The scaling of \Gh\ can be inferred by applying a renormalization group transformation 
/i = bfi' and imposing hy^rh{^) = 1 as in Sect. 13.11 The LHS becomes: 

Co(t;m,/x) = m^Co(tm^/^'";l,/i) , (49) 

whereas the RHS scales as 



\G 



H\ 



|2 



rvj 



2MhV 
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^^mH-l^/ym+.i/ym _ f^Q\ 



Combining Eqs. ( 149|) and (150|) we obtain: 

Aq-I 

ICi/l r^ m y^ . (51) 

The definitions of tlie decay constants, tlieir anomalous dimensions and resulting scaling 
coefficients are summarized in Tab. [TJ 

We would like to the draw the reader's attention to the fact that the pseudoscalar 
decay constant, as defined in Tab. [H is related through the PCAC relation as 

d-A'' = 2mP'' , ^ 2mGpa = MJ,aFpa . (52) 

The scaling is consistent with our findings from Eq. fl5T]) since: 

1 + (2 - 7,)/y„ = 2/y^ + l/vm . (53) 

Let us briefly discuss the scaling dimensions of the operators given in Tab. [H The currents 
V, V"" have vanishing anomalous dimension since they are conserved currents that are 
associated with global symmetries. The axial current A"- is only partially conserved, see 
Eq. (l52l) . It is broken by a soft term whose renormalization does not affect the divergence 
d-A"- and therefore A°- has vanishing anomalous dimension. Moreover this implies that 
mP"' is a renormalization group invariant and thus A pa = 3 — 7*. The scaling dimension 
of S was already discussed in section |2l In the case where there are no masses S"" and 
P"^ have the same renormalization constant. This is explicit to all orders in perturbation 
theory and should also hold non-perturbatively. Neglecting effects of the mass on 7* 
one concludes A^a = 3 — 7*. The flavour singlet axial vector identity is anomalous. 
The topological charge density mixes with the axial vector, which therefore does not 
renormalize multiplicatively. This is further discussed in appendix IB. 21 

4.2 Low energy theorems from Ward Identities and alike 

In the previous subsection we have inferred the scaling laws of the decay constants of the 
lowest lying states from the anomalous dimensions. Further information can be obtained 
by analyzing Wis. 

In appendix [B] we recall the derivation of two standard Wis, Eqs. flA.5p and flA.12p . 
and a low energy theorem, Eq. (1A.13P : 



(2m)^np„pj0) = -2m6ab\ 
{2myUpp{0) - Uqq{0) = -4m{qq) 

^ss{0) = --^m , (54) 
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Table 1: Scaling laws, G[F] ~ 7n'''^[^l, for the decay constants of the lowest lying states. In 
regard to the V/A decay constants and formula (fSTj) note that Gy/A ^ ^v/aFv/A in terms of 
counting scaling powers. The symbol Hm = 1 + 7* denotes the scaling dimension of the mass 
([7]). Recall that the lowest bound state scales as Mh ~ m^'y"^ p4p and the non-analytic part 
of the quark condensate is given by: (qq) ~ m^'^"'^*-''^'" (j3]). The symbol a denotes the adjoint 
flavour index, and A" are the generator normalized as tr[A°A ] = 26 . 

where UxYiQ"^) is the time ordered two-point function, c.f. (IA.6p . Information on the 
decay constants can be gained by investigating the dispersion representation of the two- 
point functions: 

Im[np,P,(.)] 



1 /" ] 

^PaP,{q') = - ds- 

TT Jcut 



+ c + dq^ + ... 



(55) 



'cut s-q^ -iO 

where we have chosen IIp^p^ as representative for definiteness. The symbols c and d denote 
subtraction constants due to UV-divergences of which only c is relevant since g^ = in 
the equations above l|. At g^ = Eq. ( 1H^ writes, 

Im[np„p,(.)] 



np„n(o) = - f ds 

TT Jcnt 



iO 



+ ci + C2m 



(56) 



where c = ci + C2m^ are the subtraction constants due to a quadratic and logarithmic 
divergence of which only ci is relevant for our discussion. 

We would like to make a comment of speculative nature. Assuming that the lowest 
state. 



np.n(o) = 5, 



ab 



y^ pa 

ML 



+ ... 



contributes to the leading scaling of the RHS of Eqs. fl54l) 

■ ■ + Cl : 



1^2 

^^-Jf pa 



Mia 



— {m) 

m 



(57) 



(5^ 



^Note c vanishes for nQQ(O) since tlie latter vanishes to all orders in perturbation because Q can be 
written as a total derivative. 
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then, using the resuhs for Mh and setting aside the issue of the subtraction constant 
Ci ~ 0{mP) for the moment, the scahng laws in Tab. [T]are reproduced for P'^ and S and 



Gl - (Gp/2m)2 _, 2(1^^ 

M2 



^%, 1 ^ ^ 1+^,; ^ (^59^ 



'p 



would follow from the pseudoscalar WI, where the decay constants Gp and Gp are defined 
as follows. 

(0|P(0)|P(p)) = Gp , (0|g(0)|P(p)) = Gp. (60) 

The result for P, if correct, is new but on the otherhand not very practical as it is the 
difference of two positive terms. Another possibility, or rather an extension of the scenario 
above, is that the hadron mass scaling is universal and that the width does not upset the 
parametric effectdj, which would result in all decay constants scaling in the same way. 
Clearly these statements above are of a speculative nature driven by the knowledge of the 
lowest lying decay constants from the Callan-Symmazik equations. It is also amusing to 
see in what way the scaling laws for those two channels turn out to be the same. One 
might wonder about the influence and the origin of the subtraction terms in the WI ( |54l) . 
They match the divergences of the quark condensate on dimensional grounds. In fact, it 
can be seen that they match exactly. Following Ref. [56] we can fix the coupling 



mq\^n) = \l^fn (61) 

(|36|) at Aqq = 3 by demanding that the deconstructed version matches 11 (0)55 in the 
region where the theory is asymptotically free. At vanishing quark mass and 0{g^) the 
time dependent scalar correlator, Il{q'^)ss, is given by: 



P3 r sds _ NcUf 

27r y s — g^ ' An 



W)ss = ^ / -^ , P3 = ^ . (62) 



We have factored out the coefficient P3, which matches the deconstructed version. The 
condensate fl55]) in the asymptotically free region, with the normalization flM]) . becomes: 

(gg) = — m — / ds-1 . (63) 

27r J 

Whence an exact matching of the divergences in the correlation function (1541) on the LHS, 
and the quark condensate on the RHS is found. Thus the UV- divergences of the quark 



^This should be true in the large TVc-limit where the width is suppressed by 1/Nc as compared to the 
mass. 



19 



condensate and the ambiguity in defining the n(g^) functions from a dispersion relation 
do match exactly. The deconstructed version (136|) of the condensate is therefore consistent 
with the WI. Note the extra factor of one half on the RHS of the non-singlet pseudoscalar 
WI is due to the normalization tr[A°A^] = 25""^ of the flavour generatorqj. 

4.3 Comparison with data 

The scaling predictions obtained above can be compared to the recent lattice data pre- 
sented in Ref. [H] for the SU(2) gauge theory with two flavors in the adjoint representa- 
tion. The dependence of the pseudoscalar decay constant on the fermion mass is reported 
in Fig. m It is clear from the plot that the non-analytic dependence of Fpa on the fermion 
mass can not be determined from current lattice data, where no curvature is visible. This 
is confirmed quantitatively by trying to fit the data to a power law dependence on the 
fermion mass. Rather than trying to determine the exponent from the fit, we keep the ex- 
ponent fixed, and fit the proportionality coefficient only. As shown in the figure, good fits 
to the data at the smaller masses can be obtained for different values of rjpp^ by adjusting 
the constant of proportionality c. 

A preliminary estimate for the mass anomalous dimension was obtained from numer- 
ical simulations using the Schrodinger functional in Ref. [16]. A value of 7* ^ 0.5 is 
compatible with the results in Ref. |46], and leads to rjcp — 1.0 and rjFp — 0.33. This is 
represented by the blue line in the figure. As discussed above, the value c can be adjusted 
to yield a good description of the data at the smaller masses. Note that we only expect the 
scaling to hold in the limit where the fermion mass goes to zero, therefore it seems natural 
to exclude the heavier points from this analysis. However further systematic uncertainty 
is introduced by the choice of the fitting range. 

The situation improves only slightly when looking at the dependence of the coupling 
Gpa on the fermion mass, which is presented in Fig. |5l A two- variable fit of the data can 
be performed in this case, and yields a scaling exponent rjCpa = 1.3(2). Note that, using 
the scaling formula in Tab. [T|, the result of the fit implies 7* = 0.30(5), which broadly 
agrees with the result of Refs. [Ml HZl US] . 

Better control over systematic errors is required in order to extract robust information 
from the scaling of the pseudoscalar decay constant. Current data can only be used to 
check the consistency with the scaling we presented above; Fig. H] clearly shows that lighter 
masses are needed in order to actually determine the scaling exponents from lattice data. 



^We have only focused on the leading quadratic divergence of the quark condensate; it would be inter- 
esting to investigate the logarithmic divergence as well in which case one could possibly learn something 



about mass correction to the deconstructed version presented in section 13.2 
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Figure 4: The pseudoscalar decay constant Fp^ as a function of tiie fermion PCAC mass am. All 
quantities are expressed in units of the lattice spacing, which is kept constant as m is varied. The 
lines represent fits of the data to power-law scaling for different values of the critical exponents 

4.4 Further remarks 

4.4.1 Scaling of the decay width 

Let us focus on a generic decay process A -^ BC, mediated by some effective Lagrangian: 

£eff= fd^xGABcA{x)B{x)C{x), (64) 



where Gabc is the ABC coupling, and A{x),B[x),C{x) are the fields creating and an- 
nihilating the states A,B,C. The fields are normalized as (0|74(0)|74(p)) = 1 etc. Let us 
now introduce three interpolating fields J a, Jb, and Jc] these are composite fields that 
have an overlap with the single particle states, e.g. like quark bilinears for the simplest 
type of mesons. At lowest order in Gabc the correlator is given by: 



{JA{x)JBiy)Jcm ~ / d'zGABc{JAix)JBiy)JciO)C,fiiz)) 



and inserting a complete set of states becomes. 



c 



GaGbG, 

MaMbMc 



VT 

17^ 



Gabc^ 



-'MA{tx-tz)-MBity~t,)~Mci-t^) 



Using the scaling laws discussed above, we find for the LHS of Eq. (l66!l : 

{JaJbJc) ~ m(^^+^«+^^)/^- , 



(65) 



(66) 



(67) 



21 



0.25 



0.20 



p 0.15 



0.10 



0.05 



O.OC 



O 12 lattice 
° 16^ lattice 



ffil 



i* 



[iii- 



).00 



0.05 0.10 

am 



Figure 5: The pseudoscalar coupling Gp^ as a function of the fermion PCAC mass am. Ah 
quantities are expressed in units of the lattice spacing, which is kept constant as m is varied. 
The line represents a fit to the data assuming the power-law behaviour described above. 

while from the RHS we obtain: 

(JaJbJc) ~ m^^^-^)/^'" m^^B-i)/ym ^{Ac-i)/s/™ m"^/^™ m^/^™ Gabc , (68) 

and therefore 

Gabc ~ m^/^™ . (69) 

The scahng of Gabc determines the scahng of the decay width for this specific channel: 



^ ' Ma 

which corresponds to the same scaling we have argued for in section [2] Eq. (fT5|) . 



(70) 



4.4.2 Heavy quarks or mass deformed conformal? 



It has been pointed out p2j that a confining theory with chiral symmetry breaking, large 
mass term and small volume could mimic a conformal theory with mass perturbation. 
Thus the question: how to distinguish a conformal theory with small mass perturbation 
from a heavy quark regime? We would like to emphasize that the scaling laws of the 
pseudoscalar decay constant are a major help in this respect. The decay constant of 
pseudoscalar meson of two heavy quarks, which we denote by bh in analogy with QCD, 
are expected to have the following scaling behaviour: 



bb 



ra 



-1/2 



G. 



bb 



m 



1/2N 



(71) 
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This follows from the two heavy-quark state behaving like a quantum mechanical bound 
state, and therefore being treated using a quark model; Eq. (171]) expresses a result found 
a long time ago by van Royen and Weisskopf [61]. This is clearly different from the 
behaviour of the decay constant of a pseudoscalar state, or any other state, in a mCGT, 
which vanishes in the limit m — )■ 0, c.f. Tab. [H 

5 Summary and conclusions 

In this paper we have accumulated a number of analytical results for mass scaling ex- 
ponents 7] of the type (|2]) for lowest state observables from Callan-Symmanzik equa- 
tions, which should help to identify the conformal window of four dimensional non- 
supersymmetric gauge theories. Possibly the clearest evidence of the conformal window 
so far comes from the vanishing of the hadron masses and decay constants in the chiral 
limit. We have identified the scaling of Fg^ ~ m~^/^ as a criterion to distinguish mCGT 
from a heavy quark regime at small volume which has been pointed out as a potential 
pitfall in identifying the conformal window [12] . 

In section 13.11 mass scaling coefficients for condensates are determined from Callan- 
Symmanzik equations of which we mention r]qq = (3 — 7*)/(l + 7*) and 77^2 = 4/(1 + 7^,). 
In section 13.21 we provide a more physical, but somewhat more heuristic, derivations of 
those results and discuss the nature of IR- and UV-regularizations within this framework. 
By generalizing the Banks-Casher relation from QCD it is shown, in section 1371] that the 
exponent of the eigenvalue density of the Dirac operator is rjqg, providing an alternative 
method for extracting 7^,. As our discussion in section [3^ shows, it is too early to draw 
conclusions on the extraction of 7^, by this method. In section H] we derived scaling laws 
for all lowest state decay constants, summarized in Tab. [H other than those affected by 
the chiral anomaly and of tensorial structure. Fitting to current data for Fp^ and Gp^ 
we find results for 7* compatible with earlier derivations of 7* ~ 0.5. Summarizing, the 
derivations indicate that lowest state observables O scale as: 

Vo = ^, (72) 

1 + 7* 

where the one-particle state, \H{p)), scaling dimension turns out to be f]\H{p)) = ~1/(1 + 
7,^). Whether or not this is true for higher states as well, remains unclear and is more of 
academic interest as higher states are difficult to assess on the lattice. The real distinction 
of higher states is presumably the decay width and the associated continuum thresholds. 
It is therefore tempting to think that in the large A'^c-hniit, where the width is supposed 
to vanish, all decay constants and masses in a specific channel scale in the same way as 
the lowest one. 
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To this end we would like to add a few comments on tentative conclusion on the scaling 
of the S'-parameter. The electroweak S'-parameter is proportional to the V-A correlator 
evaluated at zero momentum, ny_yi(0) ~ f pv-A{s)/sds, with the pion pole subtracted. A 
lattice computation of the S'-parameter for walkinig technicolour (WTC) theories would 
be phenomenologically important. Thus it is crucial to distinguish, in a parametric way, 
the regimes of WTC and mCGT|j. Assuming that the correlator is saturated by the lowest 
resonanceql^ny_A(0) ~ Fy/My—F\/M\—F'^/Ml, it is just the pion pole that serves as an 
indicator since {M^)wtc ~ 0{m) leads to lly^{0) ~ 0{m~^) and the results in Tab. [H 
imply n™5P"'"(0) ~ 0(1) . To this end let us note that, since the pion mass is the lowest 
one in the spectrum [51], a determination of Fy/F^ < 1 would imply n™5P^(0) < 0. 
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Note added: We would like to add that in our more recent work [62] we have estab- 
lished, invoking similar RG arguments, that the scaling relation for decay constants and 
masses are valid for the entire spectrum. Furthermore, using the trace anomaly and the 
Feynman-Hellmann theorem we have derived the mass scaling relations without resorting 
to RG arguments. 

A Operator product expansion in the deep Euclidian 

In Ref. [63] an OPE relation is obtained which we shall derive here in a slightly modified 
form. Taking the Fourier transform of Eq. flA.4p one arrives at the expression 



2mq^ j e''^-^{TAl{x)P'm, = {2mflip^p^{q^) + 2m{qq) . (A.l) 

Inserting a complete set of states on the LHS and using Eq. (152!) yields: 



^Of course this does not solve the problem of distiiiguishmg the WTC from TC regime per se. 
^"This is reasonably satisfied in QCD. Adding one more triplets of states P, V, A would not really alter 
the conclusions above. 
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where Q^ = —q"^, and di, 62 are subtraction constants, which can also depend on m. 
Neither of these are of relevance to us since we may simply differentiate this expression 
twice with respect to Q^ . Following Ref. [63], an expansion in one inverse power of Q^ 
yields Eq. ( I58ll . by assuming that Yip^p^^—Q"^) does not vanish as Q^ — ;■ 00 and m — )■ 0. 
By expanding in one more inverse power in Q^ one observes that the scaling of FpgMpg 
has to be larger than m^, 

2vF,s + 4Wps > 2 . (A.3) 

In the domain where 1 < 7*, Eq. (lA.Sp leads to 7* < 2, which corresponds exactly to the 
unitarity bound A^^ = 3 — 7 > 1 for a scalar field [T7]. 

B Ward Identities and a low energy theorem 

In this appendix sketch the derivation of two standard WI and one low energy theorem. 

B.l Ward identity for pseudoscalar non- flavor- singlet 

Starting with the identity, 

d^{0\TA;{x)P\0)\0) = {0\Td-A%x)P\0)\0) + 6{xo){0\[A^,{x),P\0)]\0) , (A.4) 

and integrating the equation over t/^x one arrives at the pseudoscalar Ward identity 
(WI0: 

(2m)2np^pJO) = -2m6ab{qq) , (A.5) 

where 

np.n(g) = ^ /e*'^-^(0|TP„(x)A(0)|0) . (A.6) 



Throughout the paper IlAB{q) denotes the time-ordered two-point correlator of operators 
A and B. The RHS of Eq. flA.SP is obtained by evaluating the commutator 

{0\[Qt\.o=o,d-A\0)]\0) = 2t5^'m{qq) . (A.7) 

The charge is defined as usual: Q'^I^^^q = J (PxA^{0,x). 



^We are using the fact that there are no Goldstone bosons since S\J{nf)A is exphcitly broken. 
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B.2 (Anomalous) Ward identity for pseudoscalar flavor-singlet 



The flavor singlet sector can be analysed by reviewing the Ward identities [64] used in 
discussing the ?7'-mass/U(l)yi-problem in QCD. Let us define the following fiavor singlet 
quantities: 

^(^) = Y^ Qj'^l^qjix) , A^{x) = Y^ qj7f,75qj{x) , Q{x) = -^^nfGa,pG'^^{x) (A.8) 
j=i i=i 

where Gap = '^I'^^ap-ysG'^^ . The anomaly equation is given by: 

d-A = 2mP + Q (A.9) 

The integrated anomalous Ward identity is readily obtained from Eq. ( ]A.4p and reads: 

{2mflipp{<d) + (2m)nQp(0) = -4m(gg) . (A.IO) 

By observing that 

= I f d,{0\TQiO)A'^{-x)\0) = (2m)nQp(0) + nQQ(0) , (A.ll) 

J X 

the WI f lA.lOp can be written in a more symmetric form: 

{2myUpp{0) - Uqq{0) = -4m{qq) . (A.12) 

B.3 Low energy theorem for scalar flavour-singlet 

A simple and usefel relation follows from the fact that the operator S{x) = qq{x) appears 
in the Lagrangian. Eq. ([1]) implies: 

UssiO) = --^{qq) , (A.13) 



where Ussiq"^) is the time ordered two-point function flA.6|) 



C The decay width 

In this appendix we investigate whether the large tMn behaviour of a correlation function 
can be influenced by the width. According to [65] the Euclidian time behaviour of a two- 
point function in the rest frame of the decaying particle can be written as a spectral 
integral of the type 

Cnit; g,m,fi) = ^j dEe-'''p{E) , (A.14) 
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where 

and S(-E) is the self energy and m the bare mass. We shall work in the approximation 
where E(£') does not vary appreciably around the peak E = \M\ and we neglect the 
far away singularity E = —\M\. The symbol M denotes the renormalized mass: M = 
m^ — Re(S(M)). In this case p{E) assumes the form: 

where 7 = r/2 = Im(S(ii^))/2M. Then the two-point function takes the following form 
[65]: 

Cuit; g, m, /i) = ^^Ei(7t, (M - \)t) , (A.17) 

where A is the onset of the cut, omitted in Eq. ( JA.141) . and 

Ei(a,/3)=/ — -dx>0. (A.18) 

J_I3 x^ + a^ 

If there are cancellations between the ~ e~*^* behaviour of the mass and the width then 
the following must be true: 

Ei(«, /3) *>-' «"el''l" , (A.19) 

where n a real number and 6 is a number that would have to be fined tuned. It can be 
shown that this cannot be the case. Consider 

1"^ a^e""^ f°° 

aEi{a,/3)= / — -dx < / e~''dx < 00 , (A.20) 

J ^13 x^ + a^ J_p 



but then from (1A.19P : 



cx) > aEi(a, /3) *>-' a^+^el^'l" "4°° divergent (A.21) 

one gets an immediate contradiction to the hypothesis (JA.191) . Thus we have shown that, 
in the approximations mentioned above, that a the large Euclidian time behaviour of the 
mass and the width do not conspire to cancel each other. 
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